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I. INTRODUCTION 

In this paper we answer two questions about quantum maps. 
The first question is, given any completely positive trace- 
preserving (CPT) map, how can one determine whether this 
map can be decomposed as a convex combination of unitary 
maps? In more formal terminology we ask for necessary and 
sufficient conditions such that a CPT map $ can be written as 



(1) 



where the scalars pi form a probability distribution (i.e. they 
are non-negative and add up to 1) and where the Ui are uni- 
taries. Preferably, the method fordoing so should be construc- 
tive and should as a bonus supply the pi and Ui. CPT maps 
obeying this condition are called Random Unitary maps. The 
second, and very much related, question we answer here is 
about finding a proper distance measure between a CPT map 
and the set of Random Unitary maps, and methods for calcu- 
lating it, be it numerical or (preferably) analytical ones. 

The physical motivation behind these questions is the desire 
to distinguish various error mechanisms afflicting the prepara- 
tion and processing of quantum states. For example, errors oc- 
cur in the realisation of quantum gates in quantum information 
processing. If the only error mechanism occurring is a clas- 
sical uncertainty, for instance in a phase parameter, then the 
resulting "gate" will not be described by a particular unitary, 
but rather by a mixture of such unitaries; the mathematical 
description of such a mixture is effectively a Random Unitary 
map. If on the other hand, other mechanism can occur, such as 
spontaneous decay or photon loss, then the resulting gate can 
no longer be described by such a Random Unitary map. In a 
sense, the distance between this particular map and the set of 
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Random Unitary maps determines and quantifies the presence 
of these non-classical error mechanisms. 

The paper is structured as follows. In Section|IIl we present 
notations and definitions for a number of basic concepts that 
will be needed in the rest of the paper. Section [Til] is devoted 
to the determination of whether a unital CPT map is a random 
unitary map. The techniques introduced in this Section are 
then generalised in Section[lV]to construct a genuine distance 
measure D between a map and the set of Random Unitary 
maps. A number of properties of D are subsequently derived. 
Finally, in Section [V] a connection is made to the entangle- 
ment of assistance of bipartite states. 



II. PRELIMINARIES 

It is quite obvious that any Random Unitary map should, 
apart from being trace-preserving, also be unital, meaning that 
$(1) = 1, Maps like this, trace-preserving and unital ones, 
are also called doubly stochastic. This necessary condition 
of double stochasticity implies that the input and output di- 
mensions of the map should be identical. We shall hence- 
forth assume that $ is a unital CPT map taking states on a 
d-dimensional Hilbert space 7i to states on that same space. 

At this point it is of course very tempting to check whether 
the condition of double stochasticity may even be necessary. 
In fact, for qubit CP maps (d = 2) this is the case. By defini- 
tion, the set of doubly stochastic CP maps is convex. In (jl, 
Landau and Streater proved that for d = 2 the extremal points 
of this convex set are precisely the unitary maps. This is just 
a reformulation of the statement that the set of Random Uni- 
tary qubit-maps is precisely the set of doubly stochastic CP 
qubit-maps. 

For higher dimensions this is no longer true, as was first 
shown by Tregub and by Kummerer and Maassen til |2j, |3[] . 
That is, for higher d there are extremal doubly stochastic CP 
maps that are not unitary. An example of such a map in odd 
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dimensions (taken from Ref. [1]) is the following: 

where d = 2j + 1 and the Jk are the three well-known gener- 
ators of SU (2) in its d-dimensional irreducible representation 
|0] • Consequently, the set of all convex combinations of uni- 
tary maps is a proper subset of the set of all doubly stochastic 
CP maps. 



A. Jamiotkowski Isomorphism 

To proceed, we will next exploit the Jamiolkowski isomor- 
phism between CP maps and states. For the purposes of this 
paper it is not strictly necessary to do so, but it has the benefit 
of widening the perspective. The Jamiolkowski isomorphism 
assigns to each CP map from H. m to fi ou t a (not necessarily 
normalised) state on 7i out (gi Hi n ; this state will be called here 
the state representative of the CP map, and will be denoted 
by <£> as well. The context will make clear whether the map 
is meant or its state representative. To avoid confusion, we 
will always assume that the state representative is normalised 
(i.e. it has unit trace). The density matrix representing the 
state representative is sometimes called the Choi matrix, after 
M.-D. Choi who proved that under the Jamiolkowski isomor- 
phism CP maps are mapped to positive semidefinite matrices 
[10], which is not true for any non-CP map. 

In our case, the state representative is a d x d state. To ex- 
plain how the assignment is done, we first introduce the sym- 
bol | J) for the maximally entangled state vector 

v i— l 

Then the following formula defines the Jamiolkowski isomor- 
phism: 

$ = ($®1)(|J)(J|), (3) 

where 1 stands for the identity map; so the map $ operates on 
the first copy of H, while the second copy is left untouched (is 
"operated" upon by the identity map). 

As our map is trace-preserving, its state representative is 
automatically normalised. Moreover, the reduction of the 
state representative obtained by tracing out the first copy of 
Tt yields the maximally mixed state II jd. Furthermore, uni- 
tality of the map shows itself through the fact that the second 
reduction (tracing out the second copy of H) is also the maxi- 
mally mixed state. 

Applying this game to the question under consideration 
yields that a unital CPT map is a Random Unitary map iff 
its state representative is of the form 

$ = 5^Pi(Ui®l)|/)(/|(tr/(8l). (4) 



This actually says that the state representative should be a 
mixture of maximally entangled (ME) pure states. Indeed, 
any such pure state can be obtained from the "vintage" ME 
state 1 7) by applying a local unitary on either party. The set of 
such mixtures forms, by its very definition, a convex set with 
the ME pure states as extremal points. We will denote this set 
by the symbol M. . 

Furthermore, we will denote the set of states with maxi- 
mally mixed reductions by M: 

Af:={p:Tr lP = Tr 2 p = t/d}. (5) 

Obviously, we have A4 C AT. We will see below that for 2x2 
states, M. is actually equal to TV, while this is no longer the 
case for higher-dimensional states. 

Our questions have therefore been reduced to determining 
whether a state is a mixture of ME pure states or, if that is not 
possible, how far the state is from the convex set A4. What al- 
most immediately comes to mind is the resemblance between 
these questions and the notorious questions of determining 
separability and entanglement, where the convex set under 
consideration is the convex hull of all pure product states. To 
wit, what we do in this paper is nothing but adapting mutatis 
mutandis various methods from entanglement theory to the 
problem at hand. 

In the next few paragraphs we introduce some more nota- 
tions that will simplify the rest of the presentation. 

B. Matrixiflcation 

An important linear operation one can perform on bipartite 
state vectors is the "matrixification" operation [5]. For a given 
bipartite d x d state vector | iji) , we define its matrixiflcation as 
ad x d matrix, denoted by tp, such that the following holds: 

h/>) = (^®l)|J). (6) 

Put less abstractly, the matrix %jj consists of all the vector 
entries of ip placed in a matrix frame. That is, if is 
given by l^) = Y^ij a ij\hj) tnen "0 is just the matrix a = 

In passing, we remark that the Schmidt decomposition of a 
bipartite pure state is nothing but the singular value decompo- 
sition (SVD) of the matrixified state vector, and the Schmidt 
coefficients are just the singular values. 

For every linear map $, one defines the dual map $' as that 
map for which Tr[A$(B)] = Tr[$'(A)B] for all A, B. If $ 
is CP, then so is Furthermore, if $ is trace -preserving then 
$' is unital, and if $ is unital then $' is trace-preserving. 

C. Bloch Vector Formalism 

We will also have the opportunity to use the well-known 
Bloch vector formalism for representing states. The Bloch 
vector formalism is based on the observation that the density 
matrices on H are themselves embedded in a Hilbert space, a 
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space of dimension d 2 . Because density matrices are Hermi- 
tian, this Hilbert space is real. The inner product in that space 
is given by the functional (p, a) := Tr[pa], The Bloch vector 
of a density matrix is just the vector representing that matrix 
in this Hilbert space, the entries of which of course depend on 
the choice of basis. 

The standard choice of basis elements are the Pauli matrices 
Oi (for d = 2) and generalisations thereof to higher dimen- 
sions; e.g., in d = 3 one has the Gell-Mann matrices. We use 
the notation Tj for these generalisations. They can be grouped 
into x, y and z groups H : 

T X ;kl = \k)(l\ + \l)(k\, l<k<l<d 

Ty;kl = i(\k)(l\-\l)(k\), \<k<l<d 

1 < k < d- 1. 

Here, the notation a xk stands for k entries equal to a. It is 
standard to set To = \j2/d\ . The inner products between 
these matrices are given by Tr^Tj] = 2Sij (the factor 2 is 
conventional). 

The entries of the Bloch vector in this basis are then given 
by P : = (Po,Pi, ■ ■ ■ ,Pd 2 -i) where p t := Tr[pr;]/2; con- 
versely, p = J2i Pi T i- One easily checks the relation Tr[/9er] = 
2(p,a). 

For normalised states, p = 1/ y/2d, which is a constant. 
It is therefore meaningful to employ the reduced Bloch vector 
instead, which is p := (pi,... , p^-l)- For normalised states, 
we have Tr[pcr] = 2(p,a) + l/d. 

When the state p is subjected to a unitary conjugation, p i— > 
UpW (UU^ = 1), the corresponding Bloch vector will be 
rotated according to a certain orthogonal matrix O (00 T = 
1). Let p = J2i Pin, then p' = UpW = PtUnW and 

p'j = Tr[pS-]/2 

= ^piTrpTiU^Tj}^ 

i 
i 

or p' = Op, where Oji ~ Ti[UTiWrj]/2. This defines a 
real matrix O (because Bloch vectors are real), and we need 
to show that O is orthogonal. That is easily done as follows: 

(oo T ) t3 = Y.°^o Jk 

k 

k 

= Tt ^- U fer fc Tr[[/T fe f/t Tj ]/ 2 j [/t r .]/ 2 
= Ti[U(Wt 3 U) C/ t r l ]/2 

so that, indeed, OO t — 1. Furthermore, as the unitary con- 
jugation leaves tq invariant (it is a multiple of the identity), 



O decomposes as [1] © O, where O operates on the reduced 
Bloch vector only. 

Positivity of the density matrix translates to certain condi- 
tions on the Bloch vector. It implies that the purity Tr[p 2 ] 
of states lies between l/d and 1, so that the length of the 
Bloch vector is bounded from above by 1/ \/2. Likewise, the 
maximal length of the reduced Bloch vector is -J(d— l)/2d. 
States whose reduced Bloch vector is this long are automati- 
cally pure. This means that the reduced Bloch vectors of states 
all lie in a ball of radius y/(d— l)/2d, the so-called Bloch 
ball. 

For d = 2, all points in the Bloch ball turn out to correspond 
to states. For higher d this is no longer the case, and the set 
of Bloch vectors of states is more like a "dimpled" ball, the 
conditions on positivity leading to slices being cut away of 
the original ball. It is easy to see that this has to be so: the 
vector of diagonal entries of a density matrix is a probability 
vector and thus lies on a simplex, which means that there must 
be a linear projection under which the set of all Bloch vectors 
turns into that simplex. This is impossible for a ball, unless 
that simplex is 1 -dimensional, as it is in the d = 2 case. 

III. DETERMINING WHETHER A UNITAL CPT MAP IS A 
RANDOM UNITARY MAP 

To answer the question of whether a unital CPT map is 
a Random Unitary map, we freely borrow the methods de- 
scribed in [5] and adapt it to the problem at hand. 

A. Condition for Random Unitary-ness: pure case 

The first step in answering the question is to find a quadratic 
criterion for judging whether a pure state is maximally entan- 
gled (ME). This is quite simple, and the answer is that a state 
vector is ME iff Vdip is unitary. Indeed, we noted that ME 
states are characterised as \ip) = (U (g> t)\I) for some unitary 
U, so we trivially get ip = U. 

Thus, one must have 

i>$ = l/d. (7) 

The left-hand side is easily seen to be the reduction of the state 
\ip) (ip | to the first subsystem; thus the condition can also be 
phrased as 

Tr 2 |V)(^| = l/d. (8) 

It is immediate from the purity of the state that the reduction 
to the other subsystem will automatically be 1 /d as well. For 
normalised tp, this leaves us with a system of d(d + l)/2 — 1 
quadratic equations. 

In terms of CPT maps, pure states are the state representa- 
tives of rank-1 CP maps, i.e. maps with a single element in 
their Rraus decomposition. The above discussion then leads 
to the rather obvious fact that rank-1 CP maps are Random 
Unitary maps if and only if they are unital and trace preserv- 
ing. 
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For higher-rank CP maps, this is no longer the case. Never- 
theless, the condition Eq. ([8]l will play a central role for such 
maps as well. It will turn out to be convenient to rephrase that 
condition in terms of the Bloch vector of the reduced state. 
One easily sees that it is equivalent to the condition that the 
reduced state has the zero-vector as reduced Bloch vector; this 
also holds when the state is non-normalised. Thus, the condi- 
tion is given by the system of equalities 

Trfa Tr 2 | ip) (t/>|] = 0, for 1 < i < d 2 - 1, (9) 

or 

(tp\n ® t\4>) = 0, for 1 < i < d 2 - 1. (10) 

B. Decompositions of Mixed States 

Now we must use the pure-state criterion to see whether 
a mixed state can be decomposed as a mixture of ME pure 
states. This brings us to the second step of the method: de- 
termine all possible convex decompositions of a state. The 
answer to that is well-known indeed and has been discovered 
many times. Let the eigenvalue decomposition (EVD) of the 
state p be given by p = U MJ\ where U is the unitary whose 
columns are the eigenvectors m, and A is a diagonal matrix, 
the diagonal elements of which are the eigenvalues \; these 
are non-negative and add up to 1 . Then any other convex de- 
composition of p is given by p = WMW\ where W is no 
longer unitary and need not even be square; its columns are 
the (normalised) vectors of which the decomposition is built. 
The matrix M, though, is still diagonal with non-negative di- 
agonal elements adding up to 1 ; these are the convex weights 
of the decomposition. Alternatively, these weights can be ab- 
sorbed in their corresponding vectors so that the norm squared 
of each vector then equals its weight. This gives p = ZZ\ 
with Z = WM 1 ' 2 . In this way, each convex decomposition 
of p is related to a certain "square-root" of p. 

Characterising all possible square roots of a matrix is a sim- 
ple problem in matrix analysis. Starting from the eigenvalue 
square root, p 1 ! 2 = UA 1/>2 U Jf , one generates all others by 
right-multiplying it with a right-unitary matrix T. A right- 
unitary matrix T is a non-necessarily square matrix for which 
TTt = 1 holds (but not T^T = 1, unless T is square). Thus 
one has 

Z = p 1/2 T. (11) 

The number of columns in T is the number of vectors in the 
decomposition, and is called the cardinality of the decompo- 
sition; this number should be at least as large as the rank of p, 
i.e. the number of non-zero eigenvalues. 

In quantum physics circles, the most recent and probably 
best-known incarnation of this result is the famous Hughston- 
Jozsa-Wootters (HJW) theorem [6]. The earliest occurrence is 
actually in Schrodinger's work on quantum steering [7], and it 
has been rediscovered many times by various physicists (see 
Ref. |8[] and references therein). 

We now combine the HJW theorem with the quadratic char- 
acterisation of pure ME states to obtain a method for deter- 
mining Random Unitary-ness. 



C. Criterion for Random Unitary-ness: mixed case 

Let us consider a particular convex decomposition of our 
map, or of its state representative p, that is described by the 
matrix Z. As p is a d x d state, Z has d 2 rows. The cardinality 
K of the decomposition, being the number of columns of Z, 
cannot be smaller than the rank R of p. We want to check if all 
the vectors in this decomposition are ME state vectors. Thus 
we have to apply the criterion ( TTOb to every column vector of 
Z. This yields the system of equalities 

(Z f {n <x> i)Z)jj = 0, for 1 < j < K, 1 < i < d 2 - 1. (12) 

Inserting Eq. ( fTTT i gives 

(ry/2 (Ti ^ 1)jO i/2 T)3 _ 0; (B) 

for 1 < j < K, 1 < i < d 2 — 1. For succinctness we will 
introduce the d 2 — 1 matrices A{ defined by 

At ^^(Tigljp 1 / 2 . (14) 

The problem is thus reduced to the following: 

Problem P: Find a scalar K > d 2 and a right-unitary d 2 x K 
matrix T that " off-diagonalises" the d 2 — 1 matrices Ai (of 
dimension d 2 x d 2 ) simultaneously: 

Diag(T t A,T) = 0, for 1 < i < d 2 - 1. (15) 

Necessary conditions for such a T to exist are that all Ai 
should be traceless. For the Ai of Eq. ( TBI this means that 
Tr[Ti T12 [p]] = 0, i.e. Tr 2 [p] must be proportional to the iden- 
tity, which is a condition we have already encountered. 

In general, we do not know how to solve this problem ana- 
lytically, and we have to resort to numerical methods, just as 
in I5J. This will be described below in Section[IV]on distance 
measures. 

The off-diagonalisation problem can easily be solved for a 
single matrix A. From the lemma below it follows that the 
necessary condition Tr[A] = is also sufficient. Here, K — 
d 2 suffices, and T is a genuine unitary matrix. 

Lemma 1 For a Hermitian nx n matrix A, a unitary T exists 
such that (T'AT)jj = Ofor all j if and only if Tr[ A] = 0. 

Proof. Necessity is obvious as the trace is unitarily invariant. 
To show sufficience, consider Schur's majorisation theorem 
Jll], which says that for any Hermitian matrix X, 

j'=i i=i 

and equality holding for k = n. Horn's Lemma [9] adds to 
this that for any specified set of diagonal entries and eigen- 
values obeying this majorisation relation a Hermitian matrix 
exists exhibiting those diagonal entries and eigenvalues. 

For fixed Hermitian X the mapping j ^jPO * s 
non-increasing by definition, so that the mapping k 1— > 
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Sj=i ^j(X) is concave. Now note that if Tr[X] = then Let us also put 
E"=i *jP0 = 0, while J2° 3 =i A*P0 - trivially. By the 
above concavity statement we then find that J2j=i (^0 — 
for all k. As the zero vector (0, . . . , 0) is majorised by any 
non-negative vector, Horn's Lemma then implies the existence 
of a Hermitian matrix with zero diagonal and any prescribed 
set of eigenvalues that add up to 0. 

Consider now the eigenvalues of A, which add up to by 
assumption. By the above, another matrix X must exist ex- 
hibiting the same eigenvalues, hence unitarily equivalent to 
A, and with zero diagonal. Therefore, a unitary T exists such 



that X j:j = (T^AT) 3] = for all j 



□ 



For qubit maps there are three Ai to cope with. However, 
in the light of Landau and Streater's result that the qubit Ran- 
dom Unitary maps are exactly the doubly stochastic maps, the 
condition Tr[A;] = should also be sufficient for d = 2 and 
Ai of the form < TT~4T >. Indeed, the following Theorem holds, 
which therefore supplies an alternative proof of Landau and 
Streater's Theorem. 

Theorem 1 Let p be a 2 x 2 state with reductions p\ = pi = 
1/2; let p = U AC/t l, e the eigenvalue decomposition of p; let 
W = Uh 1 / 2 . Let also <Ji be the SU(2)-Pauli matrices. Then 
the diagonal elements ofW^ai (£> 1)W are Ofor all i. 

In other words, in this case we do not even have to search for 
the matrix T. 

Proof. The first thing to note is that the diagonal elements 
of W^(ai (g> 1 ) W are if and only if the diagonal elements of 
U^(<7i ® t)U are 0. As pure states in Af are automatically in 
AA, an equivalent statement of the Theorem is that a 2-qubit 
state is in Af if and only if all its eigenvectors are. 

A 2 x 2 density matrix can be represented by a 2 x 2 block 
matrix: 



P : 



B C 
Ct D 



where every block is a 2 x 2 matrix itself. The conditions on 
p then translate to 

B + D = ^1, Tr[B] = Tr[L>] = 1/2, Tr[C*] = 0. 

If we drop the normalisation condition on p, these conditions 
relax to 

B + D = kl,Tr[B] = Tr[D],Tr[C] = 0, 

for some real number k. 

We will now show that if p satisfies these conditions, then 
its square also does. Without loss of generality we can apply 
a local unitary, so that B can be diagonalised. Let us set 



B 



a 
b 



with a, b non-negative real numbers. The conditions on p then 
demand that D is given by 



D = 



b 
a 



C 



z x 



where x, y, z are complex numbers. 
The square of p is now given by 



P 



B 2 + C0 BC + CD 

(bc + cdY c^c + d 2 



Let us now test the required conditions: 

J B 2 + D 2 + CC t +C f C = k't, 

Tr[B 2 + CC^] = Tr[D 2 + C t C], 
Tt[BC + CD] = 0. 

Straightforward calculations reveal 

B 2 + D 2 = (a 2 + 6 2 )l,CC t + C t C= (|x| 2 + |?/| 2 + 2|z| 2 )]l, 

so that the first condition is satisfied. Since CC< and C'C 
have the same trace, and Tr[i? 2 ] = Tr[L> 2 ] = a 2 + b 2 , the 
second condition is also satisfied. Finally, Tr[BC + CD] = 
(a + b)z — (a + b)z = 0, which shows that p 2 is of the required 
form. 

We can now repeat the process of squaring and find that 
any 2™-th power is of that same form. Let us now invoke 
the power method for finding the dominating eigenvector ip 
of p: for p with non-degenerate spectra, its m-th power, after 
normalisation, tends to when m tends to infinity. By 

the above, we thus find that the projector on the dominating 
eigenvector of p is also in J\f. Therefore, if we "deflate" p by 
subtracting \\ \ifj) from it and renormalise, we again obtain 
a state in Af. Continuing in this way, we thus find that every 
eigenvector of p is in Af, and hence in AA. 

For p with degenerate spectra, continuity considerations 
lead to the conclusion that one can always find vectors in its 
eigenspaces that are in AA. □ 



D. Extremal CPT and UCPT maps 



The following theorem by Choi IIHJH characterises the ex- 
tremal CPT maps. 

Theorem 2 (Choi) The CP map defined by p i— > <l>(p) = 

Sfc=i ^feP^fe is extremal within the set of CP maps with pre- 
scribed value o/$'(l) if and only if the set 

{A k Aj;k,l = l,...,R} 
of R 2 matrices is linearly independent. 

For CPT maps, the requirement is $'(1) = 1. As a sim- 
ple consequence of this Theorem we note that the extremal d- 
dimensional CPT maps have rank at most d. This is because at 
most d 2 matrices of size d x d can be linearly independent. In 
other words, the convex set of CPT maps on a d-dimensional 
Hilbert space is the convex hull of the set of CPT maps of rank 
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at most d. Numerical experiments for d up to 6 reveal that in 
d dimensions one can find d matrices Af. satisfying this con- 
dition, so that the constraint R < d is saturated. Furthermore, 
the condition is generically satisfied for a randomly generated 
set of d matrices Ak- 

Let us now proceed to the study of extremal unital CPT 
maps. The relevant generalisation of Choi's theorem is (fTj], 
Theorem 2): 

Theorem 3 (Landau-Streater) The CP map defined by p i— > 
$(p) = ^2k=l A\.pAk is extremal within the set of CP maps 
with prescribed values of^(l) and $'(!) if and only if the set 
of R 2 matrices (of size 2d x 2d) 

{A\Ai © AiA\; k,l = l,..., R}. 
is linearly independent. 

This directly implies ([1], Remark 3) that the extremal maps 
have rank not higher than \[2d. Again, numerical experiments 
for d up to 6 reveal that in d dimensions one can find [\/2dJ 
matrices Ak satisfying the condition of the Theorem, so that 
the constraint R < \[2d is saturated; moreover, the condition 
is generically satisfied for a randomly generated set of LV2dJ 
matrices Ak- 

This remains particularly true for unital CPT maps, up to 
one exception: for d = 2 one cannot find more than 1 matrix 
Ak obeying the independence condition. This is in accordance 
with the statement ([ij], Theorem 1) that for d = 2 there are 
only rank-1 extremal maps. The existence of rank-2 extremal 
unital CPT maps is prevented by the conditions for double 
stochasticity, A X A\ + A 2 A\ = 1 and A\A x + A\A 2 = 1, 
which imply that the singular value decompositions of A\ 
and A 2 must be A x = UY,{V^ and A 2 = UYs 2 V\ (with 
the same U and V\) with £ 2 + E| = 1. Therefore, the set 
{A^A; © AiA\.; k, I = 1, 2} is not independent. Indeed, two 
of its elements are equal: 

A\A 2 ®A 2 A\ = VY^V^ © C/EiSaC/t 
= A\A X © A X A\. 

In this context, the following conjecture is of relevance 
ifTTll . The conjecture is supported by numerical evidence. 

Conjecture 1 (Audenaert-Ruskai) Every d x d state p can 
be written as an equal-weight average of d states pi (not nec- 
essarily different) that are of rank at most d and have partial 
traces Tta Pi and Tr^ pi identical to those of p. 

Remark. In numerical experiments one is confronted with 
the question of how to generate random CPT maps, unital CP 
maps, and doubly stochastic CP maps. The first two questions 
are readily solved: one generates a random CP map, and then 
projects it onto the set of CPT maps or unital CP maps, respec- 
tively. Here, the two respective projections are the operations 
(performed at the level of the map's Choi matrix): 

$ i-f (Gi ® l)$(Gi ® 1), Gi = (Tra^])- 1 / 2 
$ i-f (1 ® G 2 )$(l ® G 2 ), G 2 = (Tn^])- 1 / 2 . 



Note that these projections preserve CP-ness and do not in- 
crease the rank of $. 

The question of how to generate random unital CPT maps 
is slightly harder, as one has to satisfy the two constraints of 
TP-ness and unitality at once. Fortunately, this can also be 
done using a projection method. The method, called "projec- 
tions on convex sets" (POCS), consists of an iterative scheme 
whereby the two projections G\ and G 2 are alternatingly ap- 
plied to an initial CP map. It turns out that, due to the convex- 
ity of the two sets, this process converges very quickly lfl2ll to 
a CP map in their intersection, i.e. to a CP map that is both 
CPT and unital. 



IV. DISTANCE MEASURES 

While up to this point we have been looking at conditions 
under which a unital CPT map is a Random Unitary map, we 
can modify our method slightly to calculate a kind of distance 
between a given map and the set of Random Unitary maps. In 
general, one can define distances between a point and a set as 
the minimal distance between that point and any point from 
the set. Choosing different point-to-point distance measures 
thus induces different point-to-set distance measures. The ex- 
act choice of distance measure may be guided by consider- 
ations of physical relevance or just of mathematical conve- 
nience. An example of a distance measure with a clear physi- 
cal meaning and relevance is the gate fidelity 11311 . However, 
the Hilbert-Schmidt norm distance between the state represen- 
tatives is arguably the simplest one when it comes to actually 
performing the minimisation. 

To calculate the chosen point-to-set distance, the point-to- 
point distance has to be minimised over all points of the set. 
The problem we encounter in our situation is that the set is 
defined in terms of its extremal points, of which there is an 
infinite number. The minimisation has thus to be performed 
over all possible convex combinations of an infinite number of 
points. By Caratheodory's theorem [ 14], only a finite number 
of points can have a non-zero contribution to the convex com- 
bination. If the set is embedded in a d-dimensional real space, 
the maximal number of points required is d. In the present 
case (maps on C d ), we are dealing with d 2 x d 2 PSD ma- 
trices, which are embeddable in a d 4 -dimensional real space. 
Therefore, we need at most d 4 points to make up the convex 
combination. Nevertheless, the minimisation consists of vary- 
ing d 4 — 1 real convex weights and d 2 — 1 real parameters (to 
make up an SU (d) -unitary) for each of the d 4 extremal points, 
hence of the order of d 6 parameters in total. 

In the following we take a different approach, by modifying 
the treatment from the previous Sections such that a quantity 
emerges that is more easy to calculate than induced point-to- 
set distance measures but also has an interpretation as a dis- 
tance measure. To do so, we take the vectors of diagonal el- 
ements Diag(T^iT), concatenate them into a single vector 
and then find the T that minimises a well-chosen norm of that 
vector. The minimal norm then quantifies how far the given 
map is from the set of Random Unitary maps. At this point, 
we cannot yet say that what we get in this way is a genuine 
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distance measure. What we do get already is that the given 
map is Random Unitary if and only if this minimal norm is 0. 

First of all, we have to properly choose a vector norm. We 
can, for example, choose a norm of the form 



Property 2 The value of D for pure d x d states is given by 



D P,q(P> T ) 



1/9 



33 Ji 



and its minimal value 

which is to serve as distance quantification between the given 
map and the set of Random Unitary maps. In what follows, 
we will restrict our attention to Z?2,i, which we denote by D 
without subscripts, because it has a number of desirable math- 
ematical properties. For instance, it can be redefined as the 
convex hull of a simple function on pure states (see Property 
EJbelow). 

Calculating D now requires minimisation over T, which is 
a right-unitary matrix of dimension d 2 xK with K at most d 4 ; 
again there are of the order of d e parameters, but in this case 
they are all contained in a single mathematical object, a right- 
unitary matrix, rather than in several unitaries and a number 
of convex weights. The minimisation therefore has a simpler 
mathematical structure, which leads to simplifications at the 
level of actual algorithms, but also concerning the derivation 
of its basic properties. In [5], a modified conjugated gradi- 
ent method is described for minimising functionals over the 
manifold of right-unitary matrices. This method is directly 
applicable to the problem at hand, and we have implemented 
it in Matlab ||l5ll. 

In the rest of this Section, we discuss various properties of 
D, including lower and upper bounds that are easy to calcu- 
late. 

Property 1 D(p) is invariant under "local" unitaries, that is, 
unitaries operating on input or output space separately. 

Proof. When p is subjected to local unitary rotations, p i— > 
WpW\ where W — U <B> V, the Ai matrices will transform 
according to 

Ai i-> Wp 1/2 W\n®l)Wp 1/2 W ji 
= Wp 1/2 (U' ( T l U(Z>t)p 1/2 WK 

The "outer" appearence of the unitary conjugation will of 
course be absorbed in the unitary T in (T^AiT) and plays 
no further role for determining D(p). The unitary conju- 
gation on the n corresponds to replacing n by J^k OikTk 
(1 < i < d 2 — 1), for some real orthogonal matrix O (see Sec- 
tion lll Cl >. The entries (T^AiT)jj are replaced accordingly by 
J2k Oik{T^ AkT)jj, which amounts to rotating each one of 
the vectors ((T^ AiT)jj)i (for all j). As rotations leave the 
length (^2-norm) of a vector unchanged, this shows that D(p) 
is indeed invariant under local unitaries. □ 



D(\i>)m = V2(MTr B hMVI) 2 ] - l/d) ( 16 > 
= \/2 \\Tr B \^){ip\-l/d\\ 2 . (17) 



The maximal possible value is \/2(l — 1/d), a value that is 
achieved for pure product states. The minimal possible value 
is 0, which is achieved for ME pure states. 

Note that in the context of Random Unitary maps, we are not 
directly interested in D on all possible pure states. Neverthe- 
less, it can be calculated and, moreover, it will be useful in 
what follows. 

Note also that ( ITTb applied to mixed states in Af always 
yields 0, irrespective of whether they are in AA or not. 

Proof. We already know that T will be of no influence, 
because a pure state can only be realised in one way. We can 
therefore put T = 1 . 

For p = \if})(ip\, the matrices Ai are given by 

Ai = MM(7i®lM<# 
Then the diagonal entries of Ai are given by 

(Ai)^ = MfasiM. 

Thus, using a as a shorthand notation for Ttb(\''1>)('4 > \)< 



= ||(^|(t^i)|^)) 4 ||2 

= ||(Tr[TH)*l|2 

= *\\3h 

= 2V(Tr[a 2 ] - l/d)/2 

= V2(Tr[a 2 ] - 1/d). 



□ 



Property 3 The maximal possible value of D for a state in Af 
is v /2(d - l)/d. 

An example of a state achieving this value is the state |0) (0| ® 
t/d, which is the state representative of the CP map repre- 
senting photon loss as well as spontaneous decay of excited 
atoms: $(•)= Tr[-]|0>(0|. 

Proof. For convenience, we consider the state \d)(d\ <g> t/d 
instead. The corresponding Ai are given by 

Ai = (|d)(d|®l)(ri®l)(|d)(d|®l)/d 
d 

The only Tj with non-vanishing (d, d)-entry is r z; d-i- Its 
(d, d)-entry is given by —yj2(d — l)/d. Thus, we get 

D(p,T) = y/2(d-l)/d £|(Tt(|d)(d|® 1/4)2^1 



= y / 2(d- l)/d Tr[T\\d)(d\ <g> t/d)T] 
= y / 2(d-l)/d Tr[|d)(d| ® 1/d] 
= V2(d-l)/d. 
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As this value is independent of T, this is also the value of 
D{p). 

By Property 12 this is the maximally achievable value of D 
throughout, and what we have just shown is that this value is 
achievable for states in TV. □ 

Property 4 The function D is the convex hull of its restriction 
to pure states. That is: 

D(p) = min/^ ftj D(|^)<^|) : 

V J 

As a direct consequence, D is a convex function. 

Proof. We proceed in much the same way as we did in Section 
MI CI where we obtained a criterion for Random Unitary-ness 
in the mixed case. Again we need to consider all possible 
ensembles realising the state p, which we can do by varying 
over all right-unitaries T in the expression Z = p x l 2 T (where 
p = ZZ^). Recall that the j-th column of Z is then Zj := 
^Jpj\^j). By this variation over T we minimise the quantity 

3 

= Y.p^A{n®m j )) i \\ 2 

3 
3 

= Y,\W iA i T )»)ih 

3 

= D(p,T). 

The minimisation over T appearing in the convex hull con- 
struction thus, indeed, yields D (p). □ 

Remark. Since the convex hull construction has a dual, we 
can find an expression of D involving a maximisation. The 
convex hull conv(/) of a function / can be expressed as the 
Legendre transform of the Legendre transform of / [ 14]. That 
is, conv(/) = /**, where the Legendre transform /* of / is, 
in the quantum context, defined as 

r(X)=maxTr[Xp]-/(p), 
p 

where the maximisation is over all states p, and the argument 
X is a Hermitian operator. In particular, we have 

D(p) = maxTr[Xp] -D*(X) 

D*(X) = max{ip\X\ip) - D(t()), 

where the last maximisation is over pure states because D is 
the convex hull of D restricted to the pure states. 



Property 5 A lower bound on D{p) can be given in terms of 
the 2-norm distance of p to M: 

D 2 (p,M) :=mm{\\p-a\\ 2 : a € M} . (19) 

(7 

Namely, 

Proof. An upper bound on any norm distance of p to M. is 
obtained by restricting a to the extremal points of M., which 
are the ME pure states. Recall that these can be parameterised 
as | (f>) = (U ® 1)| I) for unitary U. 

Le us confine attention to pure p = first. Thus we 

are looking now at the minimisation 

min| || \^)(^\-(U(g)i)\I)(I\(U(g)t^\\ \. 

Let us now consider the 2-norm distance D2 ■ Taking its square 
yields 

D 2 {ijj) 2 = min2-2|(V>|(tf<8> 1)|/)| 2 
= 2 - 2(max | Tr[ftu] \/Vd) 2 

= 2-2w;/d. 

Note that normalisation of tp amounts to ||^||2 = 1. On the 
other hand, 

D{^) 2 = 2\\Tr B \iP)^\-t/df 2 
= 2||^t_ 1 / d ||a 

= 2{U\\\-l/d). 

A simple application of Holder's inequality yields the fol- 
lowing: for X such that ||X|| 2 = 1, ||X|| 4 > H^H^ 172 - 
Equality is achieved for X — td/\fd- 

Thus we obtain, for pure states, 

D = ^(Hi-i/d) 

> ^2(iHir 2 -iA0 

= ^(4/d)(l/(2-Di)-l/2). 

It is easily checked that the RHS is a convex, increasing func- 
tion of D 2 . Since D 2 (p) is by definition a convex function of 
p, it follows that the RHS is a convex function of p too. Now, 
D is the convex hull of the restriction of D to pure states, 
which means that the largest convex function that coincides 
with D on pure states is D itself. Therefore, the above in- 
equality on pure states readily extends to mixed states. □ 

For small D 2 , this lower bound simplifies to D(p) > 
D 2 (p) / \fd. An important consequence is that D has the de- 
sirable property of being approximately linear for states very 
close to AL just like the distance measures Di and D 2 . This 
tells us that D can indeed be interpreted as a distance (and not, 
for example, a power of a distance). 
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Property 6 A lower bound on D(p) is given by D(p) > 
max(V2|| Tr A p - l/d\\ 2 , y/2\\ Tr B p - l/d\\ 2 ). 

Proof. This is a simple consequence of PropertyH] By defini- 
tion, the convex hull of a function / that is defined on the pure 
states is the largest convex function coinciding with / on the 
pure states. The two functions appearing in the maximum are 
both convex functions that coincide with D on the pure states, 
and must therefore be smaller than or equal to D. □ 



For 2-qubit states that have rank 2, and hence for 2-element 
qubit maps, it turns out that equality holds so that we get a 
closed-form analytic expression for D: 

Property 7 For 2-qubit states p of rank 2, 

D(p) = V2max(|| Tr A p - l/2\\ 2 , || Tr B p - l/2\\ 2 ). 

While we have not been able to prove this yet, numerical 
experiments indicate that the expression also holds for 2-qubit 
states of arbitrary rank. Furthermore, this would be a simple 
consequence of the Audenaert-Ruskai conjecture for d = 2, 
combined with convexity of D and the statement for rank-2 
states. 

Proof. Let p be a rank 2 state, with eigenvalue decompo- 
sition p = p\tp)(ip\ + (1 — p)\(p)((j)\. To cover all its re- 
alising ensembles of cardinality 2, we have to consider all 

e lx cos9 sin 9 \ „ 7 , 

■a -iv a ■ We tnen § et 

— sin 9 e x cos 9 J ° 

P — \ipi) ("01 1 + iV'z) ("02 1 1 with the non-normalised states 



2x2 unitaries T = 



V>2 



e lx cos 
sin 



9\/p'4 } — sin 9y/ 1 — p (f>, 



sinOy/pijj + e %x cosflyl —p<j>. 
We then have to calculate 



£> = min\/2 V ||T(Tr A |^)(^ 



2 • 



where, in order to simplify notations, we have introduced the 
shorthand T(p) := p — Tr[p]t/d for the traceless part of a 
(non-normalised) d-dimensional state. 

Since D is invariant under local unitaries, we can take ip in 
Schmidt-diagonal form and put 

ip = (cos a, 0, 0, sina) T , 



with < a < 7r/4. As . 
form 



' is orthogonal to ip, it must be of the 



4> — (\A*sina, \f\ — rx, Vl — ry, —*Jrcosa) T , 

with < r < 1, x = sin/3, and y = cos/Je" 7 , for < < 
tt/2. One checks that |j Tr A p- 1/2|| 2 > \\Tr B p - l/2|| 2 if 
and only if |x| < \y\, or sin (3 < cos (3, or < j3 < ir/4. 

We will show below that, if < (3 < tt/4, then there exists 
a 9 and x such that T(Tr B \il>i)(i/>i\) - sT(Tr B |V 2 )(^ 2 |), 
for some s > 0. This implies that, for those 9, x, 

J2 ||r(Tr B |^)(^|)|| 2 

i=l,2 



r(Tr B |Vi)(Vi 



= ||T(Tr BP )|| 2 , 

so that D(p) < V2 \\T(Tr B p) || 2 . If, on the other hand, 
M > \vV then the same equality can be made to satisfy with 
Tr B replaced by Tr A , yielding D(p) < V% 1 1 T(Tr A p) \ | 2 . 
Now define 

K = pTCfrslVWI) 

G = (l-p)T(Tr s |0)(0|) 

H = VpU^P) (e^Tr B (|^)(0|)+h.c.) 

then 

T(Tr B IV'iXV'il) = cos 2 9 K + sin 2 9 G - sin 9 cos 9 H 
T(Tr B 1^2) (-02 1) = sin 2 6» K + cos 2 (9 G + sin 9 cos (9 H. 



We 



have 

cos 



to find 9, x such that 



sin 9cos9H = s (sin 2 (9 K + cos 2 6> G + sin 9 cos 6» H ) for 
some ,s > 0. Putting (s — 1 )/ (s + 1 ) = cos 7, this is equivalent 
to finding 9,x,7 such that 

cos{29)(K -G)= cos-y(K + G) + sin(20)#. (21) 

Now if, G and H are 2x2 traceless Hermitian matrices, 
and are determined by three real parameters: the entries of 
their respective Bloch vectors. The previous equation can thus 
be put in vector form. In order for it to have a solution for 
9, the three Bloch vectors of K, G and H must be linearly 
dependent, with real proportionality constants. 



Taking into account the special forms of ip and 4>, we have 

cos 2 a -1/2 



K = p 
G = (1-p) 



sin 2 a -1/2 

r sin 2 a + (1 — r)|a;| 2 — 1/2 \/r(l — r) (sin ay - 



Y / r(l — r) (sin ay — cosax) rcos 2 a + (1 — r)\y\ 2 — 1/2 

^ _ i—jz r / -y/r sin(2a) cos x \/l — r(cosae Ix y + siiiae^^cc) \ 

^ VI — r{cosae~ %x y + sinae lx x) — v 7 ^ sin(2a) cos x J 
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Expressing the conditions for linear dependence is the first 
step to solving our problem. Since the only non-zero com- 
ponent of K is its z-component, the three Bloch vectors are 
linearly dependent with real proportionality constants if and 
only if the (l,2)-entries of G and H (corresponding to their x 
and y Bloch vector entries) have the same argument (modulo 
7r). That is: 

Z (sin ay— cos ax) — Z(cosae lx y+sinae~ lx x) (mod it). 

The only unknown here is x> an d it turns out that there always 
is a solution. Therefore, the condition of linear dependence 
fixes x- To wit, the solution is 

e' ix = (sin(2a) — sin(2/3) cos r\ — i sin(2/3) sin 77) /Q, 
Q = y / sin 2 (2a) + sin 2 (2/3) - 2 sin(2a) sin(2/3) cos 77. 

With this choice of x> we thus nave a solution to the equa- 
tion b(K — G) = (K + G) + aH in a and b. Now we must 
make sure that this solution satisfies a = sin(20)/ C0S7, b = 
cos(2#) / cos 7 for some 7, 9. This is so provided a 2 + b 2 > 1. 
Considering only the (1, 1) -entries and the imaginary part of 
the (1, 2)-entries of K, G, H, and inserting the solution of x, 
we get a 2 x 2 system of equalities 

-H u a+(K U -G u )b = K n + G n , 
-SHi 2 a - SGub = 3Gi 2 , 

where 

K n = p(cos 2 a- 1/2), 

Gn = (l-p)(rsin 2 a+(l-r)sin 2 /3-l/2), 

/—, r r- ■ ,„ , sin(2/3) cos 77 — sin(2ci!) 

H u = y/p(l-p)y/tsm(2a)—^ ^ — , 

5Gi2 = — (1 — p)\fr sin a cos /3 sin^Vl — r, 

c^tt r~r, rCos(2a) + cos(2^) . A 

^sHi2 = V p(l — p) — cos p sin 77V 1 — r. 

We have been able to show that the solution (a, b) of this sys- 
tem, with these rather formidable expressions for the coeffi- 
cients, indeed satisfies a 2 + b 2 < 1 for < f3 < it/ A, and 
a 2 + b 2 > 1 for 7r/4 < j3 < tt/2. This was done using a com- 
puter algebra system, and we refer the interested reader to the 
supplementary material lfl5ll . This proves the Proposition for 
rank 2 states. □ 



V. ENTANGLEMENT OF ASSISTANCE 

In this final Section, we briefly touch upon another ap- 
proach to characterise Random Unitary maps which is still 
similar in spirit to the one discussed before. Rather than using 
quadratic relations to characterise whether a pure state is max- 
imally entangled, one can look at the entropy of entanglement 
of the state. Applying a similar procedure as above then yields 
the so-called entanglement of assistance. The entanglement of 



assistance, Ea, of a state lfl6ll is in some sense the converse 
of the entanglement of formation, Ep, It is defined as 

E A (p) = max { VftBfi/),) : VV^K^I = p}, 
WO V i 

(22) 

whereas replacing the maximisation by a minimisation yields 
the entanglement of formation. In other words, while Ep is 
the convex hull of the pure state entanglement functional E, 
Ea is its concave hull. 

The connection between this quantity and the problem of 
determining Random Unitaryness of a map is quite clear: a d- 
dimensional CP map is a Random Unitary map if and only if 
its state representative is a convex combination of maximally 
entangled pure states, if and only if that state representative 
has the maximal possible E A of log d. 

Fo qubit maps one has the alternative of using the concur- 
rence of assistance, Ca, defined in [17] as the concave hull of 
the pure state concurrence. Thus, a 2-dimensional CP map is 
a Random Unitary map if and only if its state representative 
has the maximal possible Ca of 2. A closed-form expression 
for Ca is given by Inll 

CA{p) = \\p 1/2 {a y ®a y )p 1 / 2 \\ 1 . (23) 



VI. CONCLUSIONS 

The set of completely positive maps can be distinguished 
by their decomposability into a convex combination of uni- 
tary conjugations. We have shown that one can find necessary 
and sufficient conditions for determining whether a given CPT 
map belongs to the set of Random Unitary maps. The criterion 
requires the collective "off-diagonalization" of d 2 — 1 matri- 
ces of dimension d 2 x d 2 that are built from the d 2 — 1 basis 
elements of the Hilbert space into which the density matrices 
on Tt are embedded. 

Based on this criterion, we have defined a proper distance 
measure to the set of Random Unitary matrices which can be 
used to quantify the extent to which non-classical error mech- 
anisms have influenced the evolution of a quantum system. 
In this way, it would be possible to point towards the dom- 
inating error source in a specific physical realization of the 
CPT map. This means that tomographic reconstructions of 
physical processes can provide a host of valuable information 
about the process itself, of which the amount of non-classical 
error mechanisms is but one. We believe that a thorough in- 
vestigation of tomographic process reconstruction can reveal 
a plethora of information about the underlying physical mech- 
anisms that led to the realization of the CPT map under inves- 
tigation, and which has hitherto not been fully appreciated. 
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